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Abstract 

The possibility to produce superheavy elements in the course of low-yield nu- 
clear explosions is analyzed within a simple kinetic model which includes neutron 
capture, 7-emission, fission and particle evaporation from excited nuclei. We have 
calculated average numbers of absorbed neutrons as well as mass distributions of U 
and Cm nuclei exposed to an impulsive neutron flux as functions of its duration. It 
is demonstrated that detectable amounts of heavy nuclei absorbing from 20 to 60 
neutrons may be produced in this process. According to an optimistic scenario, af- 
ter multiple /3-decay such nuclei may reach the long-living elements of the predicted 
"island of stability" . 

PACS: 26.30.-k, 28.70.+y, 27.90.+b, 28.20.Fc 

1. Introduction 

Synthesis of new elements is one of the main goals of nuclear physics during the last 
century. First elements heavier than uranium were synthesized in nuclear reactors via 
neutron-capture reactions followed by /3-decay. These reactions were efficient enough to 
produce elements up to Z=100 (Fm). Nuclei with larger Z are produced with accelerator- 
based experiments, most often via the fusion reaction involving a-particles and heavy ions. 
In heavy-ion induced fusion reactions new nuclei up to Z=118 were synthesized during 
the last decades [TJ [2j [3j. However, the nuclei produced in the cold fusion reactions [TJ [2] 
are situated just along the proton drip-line, thus being very neutron-deficient and short- 
lived. Fusion of actinides with 48 Ca leads to more neutron- rich superheavy nuclei with 
longer half-lives [3]. However, they are still far from the predicted "island of stability", 
where the expected life times of even neutron- richer isotopes should be much longer |Q. 
Indeed, life times of newly synthesized heavy elements with Z = 108-118, increase with 
increasing neutron content [3J. This is consistent with theoretical predictions of the island 
of stability at the neutron number N^184 and proton numbers around Z=114, Z=120, 
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and eventually also Z=126. Even more neutron-rich nuclei could be produced in multi- 
nucleon transfer processes in low-energy collisions of heavy actinide nuclei like U+Cm [5]. 
However, the cross sections of such processes are rather low. 

In this paper we discuss another possibility for synthesizing heavy and superheavy 
elements which should be more appropriate for obtaining neutron-rich nuclei up to their 
drip-lines. Namely, we propose to use intensive neutron fluxes generated by nuclear ex- 
plosions. Actually, this method was already partially employed alongside with nuclear 
reactors in the 60-s [BJ [7j. However, because of some technical and political constraints 
it was abandoned later. Nevertheless, this method provides the highest neutron densities 
which are not possible to reach with other terrestrial techniques. 

Our general ideas are illustrated in Fig. 1. It shows the upper end of the nuclear chart 
with the hypothetical "island of stability". The typical fusion reactions like 238 U + 48 Ca 
lead to proton-rich compound nuclei which lie above the /3-st ability line. On the other 
hand, the multiple neutron-capture reactions in nuclear explosions may bring the 238 U 
nuclei to the neutron drip-line. Then, after the neutron flux ceased, these nuclei undergo 
multiple /3-decay and come close to the island of stability from the neutron-rich side. If a 
single explosion is insufficient, then, in principle, a properly delayed second explosion can 
do the job. 

Several underground explosions carried out by scientists at the Livermore and Los 
Alamos Laboratories were dedicated to synthesizing new elements [Tl [8] . They include 
Par (1964), Barbel (1964), Tweed (1965), Vulcan (1966), and Cyclamen (1966). The 
elements up to 257 Fm were found in explosion debris by radiochemical methods. From 
published descriptions of these experiments one can conclude that little special effort was 
made to increase the neutron flux and irradiation time for target nuclei as compared with 
standard explosions. Also it is unlikely that the debris material, which was extracted 
for the analysis, was located near the central zone characterized by the most intensive 
neutron exposure [9] . Although the experimental equipment available at that time was not 
very sophisticated, it was sufficient to demonstrate the high efficiency of this method by 
observing the absorption of up to 20 neutrons by the U target. As discussed in refs. [HI EE] 
more heavy elements might have been produced, but they might decay via spontaneous 
fission avoiding detection with radiochemical methods. These results were discussed in the 
review article by Glenn Seaborg [6] with the conclusion that even higher-mass nuclei may 
be produced in explosions, if the neutron exposures were increased, and special targets 
were used. 

As will be demonstrated below, the multiple neutron absorption by heavy nuclei in 
the course of a nuclear explosion can lead to very high mass numbers close to the neu- 
tron drip-line. Not only a few but macroscopic amounts of such neutron-rich heavy and 
super-heavy nuclei may be produced in this way. Studying these processes may lead to 
better understanding of nucleosynthesis in r-processes associated with supernova explo- 
sions. We believe that some open questions of cosmic nucleosynthesis can be tested in 
nuclear explosion experiments by using special target materials and detection techniques. 
In our estimates below we use only open publications and Internet resources regarding 
the characteristics of explosive devices. 
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2. Main physical processes contributing to explosive nucleosynthesis. 



For our calculations we adopt a very simple model of the active zone associated with 
a low-yield nuclear explosion [TT]. We assume that it can be represented as a medium 
consisting of fissioning nuclei, fission products, light charged particles, neutrons, and 
photons. This medium is characterized by a certain temperature T which has typical 
values of 5 -j- 10 keV at the pre-expansion stage. Then we consider the time evolution 
of mass and charge distributions of nuclear species by taking into account all possible 
reactions with surrounding particles and subsequent decays. Elastic neutron collisions 
with nuclei, which do not change mass and charge distributions, are not considered here. 
It is assumed that after each inelastic interaction the nucleus relaxes to a new equilibrium 
state (compound nucleus). In particular, the (n, 7) reaction is considered as a two-step 
process, including the neutron absorption at the first step and the photon emission at the 
second step. Specifically, in our analysis we include the following processes: 

Absorption of particles by a nucleus. The neutron capture by the nucleus with mass 
number A is characterized by the width 



Here v re i is relative nA velocity which is close to the thermal neutron velocity v n = 



^J(3T/m N ) (in units c, m N = 939 MeV, mass of nucleon), p n is the density of neutrons. 
For the neutron capture cross section a n we take the inelastic part of the nA cross section, 
i.e. a n = (j t ot _ Cei- The corresponding evaluated data were taken from the web site [T2"] . 
Our analysis shows that in the neutron energy range between 2 MeV and 10 keV these 
cross sections may change by several times. Also, they show a considerable odd-even effect 
as functions of mass. For U isotopes the changes are between 1.5 barn for odd isotopes 
and 0.3 barn for even isotopes. For Cm isotopes the corresponding changes are between 
3 barn and 0.5 barn. For our rough estimates below we take o~ n =l barn for U and 2 barn 
for Cm, independent of neutron energy and mass number of the nucleus. 

At relatively low temperatures considered here the cross-sections for capturing free 
protons and light charged clusters like a-particles are negligible. However, in thermonu- 
clear explosions the 14 MeV neutrons from the fusion reaction (d + t — > 4 He+n) may 
transfer energy to d and t nuclei, sufficient to induce (d, 7), (t, 7) and other reactions [TP"] , 
which we do not consider here. 

The capture of a photon by the nucleus is characterized by the width 



Here E 1 = T(7r 4 /(30 • 1.202)) « 2.7T is the average energy of photons in matter, v 1 ~ c 
is the in-medium velocity of photons, and the density of photons is given by 



= h(v rel a n )p,, 



(1) 




r7 = h^a^E^p, 



(2) 
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(3) 
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Emission of particles from the nucleus. The width for evaporation of a particle j 
(J = n,p, d, t, 3 He, a) from the compound nucleus (A, Z) is given by: 

r, = r Xz ~ Bj ^m p " AE K B r E) ^E. ( 4) 

3 JO A ' PAZ(E* AZ ) 1 ' 

Here p AZ and p A > z > are the level densities of the initial (A, Z) and final (A' , Z') compound 
nuclei, E* AZ is the excitation energy of the initial nucleus, E is the kinetic energy of an 
emitted particle, gj = (2s 3 - + l) is its spin degeneracy factor, fij and Bj stand, respectively, 
for the reduced mass and separation energy of the ejectile. The cross section Cj(E) of 
the inverse reaction (A , Z ) + j = (A, Z) is calculated using the optical model with 
corresponding nucleus-nucleus potential [13] . The evaporation process was simulated by 
the Monte Carlo method and the conservation of energy and momentum was strictly 
controlled at each emission step. 

In the particular case of photon evaporation from an excited nucleus the corresponding 
width is given by [H] 

Jo Tc 2 c 2 h p AZ {E* AZ ) 

where E is the photon energy. The integration is performed numerically. Within the 
dipole approximation the photo-absorption cross section is expressed as: 

ctqE^T^^ 

Here the empirical parameters of the giant dipole resonance have the values o"o = 2.5A 
mb, = 0.3Er, and Er = 40.3/A 1 / 5 MeV. In the situation which we consider here the 
photon energies are in the range from about 4 MeV to 100 keV, i.e. considerably smaller 
than the resonance energy, Er w 13.4 MeV for A = 250. Thus, the photo-absorption 
cross section is rather small, cx 7 < 0.01<7o. 

Fission is an important de-excitation channel for heavy nuclei (A > 200), which can 
compete with particle emission. It is also included in the Monte-Carlo simulations. Fol- 
lowing the Bohr- Wheeler statistical approach we assume that the partial width for the 
fission of compound nucleus is proportional to the level density p sp at the saddle point of 
the fission barrier [15] : 



1 r E *Az" B f 

r / = 7T- , F » v / Ps P (E AZ -B f - E)dE, (7) 



2tt Paz (E a 



AZ) 



where Bf is the height of the fission barrier which is determined according to the Myers- 
Swiatecki prescription [16J. For p sp we have used approximations obtained in ref. [H] 
from the extensive analysis of nuclear fissility and T n /T f branching ratios. 

Beta-decay. There are many approaches to such a complicated problem as /3-decay. 
In our calculations we are mostly dealing with nuclei far from the /3 stability region. 
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Their masses and /3-decay rates have significant uncertainties. On the other hand, it is 
well known that /3-decay is rather slow in comparison with other nuclear processes, it 
has characteristic times of milliseconds and more. A maximum effect of /3-decay was 
estimated by taking the rates for 'allowed' transitions from the Sargent diagram [T7]. In 
particular, we have used the following interpolation of the decay rate A (in 1/sec): 

log{\) = Alog{E max ) - 3, (8) 

where E max (in MeV) is the maximum energy of emitted electrons, assuming that the 
daughter nucleus is in its ground state. Then the corresponding decay width is Tp = hX. 
For nuclei in the vicinity of the drip-line this formula predicts life times of order of 
milliseconds [18]. 

Spontaneous fission and a-decay. Spontaneous fission is characterized by even longer 
times than /3-decay. It will certainly affect the yield of superheavy elements before extrac- 
tion of a sample for analysis. However, it has practically no influence on the accumulation 
of neutrons by the target nucleus. To take into account this process we have employed 
a simple parametrization of the fission half-lifes proposed in ref. [8j. Typical half-lives 
for alpha-decay of these neutron-rich nuclei are usually longer than years and this decay 
mode can be completely ignored in the calculations. 

3. Physical conditions in the active zone 

Neutron density: For our further calculations we should know the density of neutrons 
generated in the active zone of an explosion. Let us assume a 100 kt explosion which 
corresponds to fissioning of 1.45-10 25 nuclei of 235 U (or 239 Pu) with total mass of 5.6 kg. 
At normal density this amount of U or Pu has the volume of about 300 cm 3 . Assuming 
that about 1.3 new neutrons are produced in average per fission, we can estimate the 
total number of free neutrons to be 2-10 25 . This gives the maximum neutron density of 
about 6-10 22 cm -3 . However, due to a finite efficiency, the amount of the fissile material 
and accordingly its volume should be larger. On the other hand, with the implosion 
assembly method, the fissile material is initially compressed by several times, that will 
increase accordingly the neutron density. According to ref. [TT] , the neutron multiplication 
process should be over within the time interval of about 1.5- 10~ 7 sec after beginning of 
the chain reaction. To this time the density of free neutrons increases exponentially and 
reaches the maximum values estimated above. Many of these fission neutrons may leave 
the core. However, some of them might be reflected back with an appropriate reflector to 
keep the neutron density in the core higher. These neutrons have a chance to slow down 
to thermal energies corresponding to temperatures T ~ 10 keV at the early expansion 
stage. Taking into account all these uncertainties as well as strong time dependence of 
the neutron yield, in our estimates below we use several values for the neutron density, 
from 4 • 10 22 cm -3 to 4 • 10 20 cm" 3 . 

One can also consider a thermonuclear reaction as a source of additional fast neutrons 
in combined fission-fusion set-ups. We have found that the neutron densities at the end of 
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the fusion stage are rather similar to the values estimated above. However, the subsequent 
evolution of the neutron density will be different in these two cases, because of a very 
different composition of the active zone. 

Characteristic time scales: There are several important time scales which influence 
strongly the efficiency of the neutron capture reaction. Obviously, the goal is to maximize 
the time r of the target exposure to the strongest neutron flux. The shortest time scale 
is the neutron multiplication time r m associated with the chain reaction. This time is 
determined by the properties of the fissile material and is typically in the range of 0.1 
fisec, as mentioned above. Such short times are obtained when all available neutrons are 
used for fission reactions, to maximize the yield of the explosion. However, this time will 
increase significantly if the active zone would contain more nuclei absorbing neutrons, like 
238 U (reduced enrichment), bringing the system closer to the critical regime. 

The disassembly stage is defined by a condition that the chain reaction is ceased. 
This happens when most of the fission energy is deposited in the active zone and the 
matter is transformed into a radiation-dominated hot and dense plasma. Due to the 
huge internal pressure it would rapidly expand into the vacuum, if nothing is done to 
confine it. However, by introducing a heavy tamper around the active zone, one can 
delay the fast expansion of plasma to a few microseconds. This time is needed for the 
shock wave to propagate through the thick layer of heavy tamper's material, usually 238 U. 
The hydrodynamic time can be estimated as a time required for the rarefaction wave to 
reach the center of the active zone, Th ~ R/c s < 10/isec, where R is the core radius and c s 
is the sound velocity in the core. The tamper can also serve as a reflector and absorber 
for the neutrons. To explore the whole parameter space, in our calculations we consider 
the exposure times r from 0.1 to 3 //sec. 

Another important time scale is the neutron thermalization time Tth in the medium. If 
the medium consists mostly of heavy nuclei like U, the neutron energy loss is dominated 
by inelastic collisions with these nuclei, followed by the photon emission from the excited 
nuclei. For slowing down the fission neutrons to thermal velocities v n (T) several (~ 10) 
such collisions are needed. Therefore, we can estimate thermalization time as Tth ~ 
10X/v n , where A is the neutron mean free path (few cm). Then one obtains Tth ~ 1/isec, 
which is usually shorter than the hydrodynamic time. 

4. Synthesis of nuclei by multiple neutron capture 

Evolution of average mass number: We have performed calculations for the 238 U and 
248 Cm target nuclei embedded in the active zone of the explosion for several fixed values of 
neutron density p n and temperature T. The 238 U nuclei are usually present in the explosive 
devices in the fissile material or in a tamper, but 248 Cm nuclei can be additionally installed 
in the coreJl] Instead of solving a complicated set of rate equations for the ensemble of 
different nuclear spices, we first consider the fate of one individual target nucleus in this 
environment. We start at t=0 with a nucleus A in the ground state and consider all 

lr The Cm isotopes are produced in significant quantities in nuclear reactors, and they have a sufficient 
lifetime to be used as a target for the neutron irradiation. 
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processes which can change its mass number or excitation energy. At each time step we 
evaluate the probabilities of these processes, pj = Tj/ (J2k^k), in accordance with the 
rates Tj discussed in the previous section. Then we choose the reaction by the Monte 
Carlo method and change accordingly the state of the nucleus. If process j is chosen, 
the clock is changed by time interval At = h/Tj. This method allows us to follow the 
stochastic evolution of an "average" nucleus. 

The average mass numbers of U and Cm isotopes are shown in Fig. 2 as functions 
of time. The dotted lines show the calculations for fast fission neutrons (mean kinetic 
energy of 1.5 MeV) assuming no moderation at all. The other lines show the results 
for thermalized neutrons at temperature T=10 keV and several neutron densities. Due 
to a higher flux, the fast neutrons dominate at early times, t < 0.3 /isec. However 
during this time interval the target nuclei can absorb only a few neutrons. The efficient 
nucleosynthesis may proceed only at longer times, when the neutrons have already slowed 
down. Then, due to the lower excitation energy, the probability of channels with re- 
emission of neutrons from the compound nucleus is strongly reduced as compared with 
reactions induced by MeV-neutrons. This is why with inceasing A the absorption rate 
becomes higher for keV-neutrons. The calculations show that a considerable part of 
curium (more than 90%) undergoes fission in reactions with fast neutrons. However, in 
the case of thermalized neutrons this channel is less probable and the majority of nuclei 
(nearly 95%) survive. One can see from the figure that U and Cm nuclei can in average 
absorb 10-15 neutrons during the time interval of 1 /xsec, if the neutron density remains 
high, 4-10 22 cm" 3 or higher. Since the neutron capture process has a stochastic character, 
even much larger numbers of neutrons could be accumulated with certain probabilities 
(see the next subsection). 

The main mechanism of increasing the nuclear mass number works as following. A 
neutron captured by a nucleus with mass A brings a few MeV excitation to the new 
nucleus A+l. This excitation energy is only slightly above the neuteron separation energy, 
which is expected to decrease with the neutron excess. It is typically below the fission 
barrier, especially in neutron-rich isotopes. The probability to return a neutron in the 
evaporation process is also low because of the small phase space volume of the final state. 
In this situation 7-emission is the main decay mechanism for reducing the excitation 
energy. After the first photon carries away some energy, the next 7-emissions become 
even more probable in comparison with other decay channels. This continues until the 
excitation is decreased nearly to zero, and then the relative probability to capture a 
new neutron becomes high again. With fixed Z and increasing A the fission probability 
decreases, since the fissility parameter Z 2 /A becomes lower. In the time interval < 1 msec 
the /9-decay still has a very low probability in order to compete with other processes. In 
the long run the survival of produced superheavy nuclei will depend on whether they can 
avoid spontaneous fission and a-decay in order to reach the island of stability via multiple 
/9-decays. 

As our calculations show, the outcome of explosive nucleosynthesis depends crucially 
on the masses and level densities of nuclei far away from the /^-stability line. There are no 
experimental data in this region, and theoretical models often give controversial results. 
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Our previous calculations were performed with Myers-Swiatecki (MS) mass formula |16j . 
The predictions based on the Cameron mass formula [19] are shown in Fig. 3 for U and 
Cm targets too. In principal, if the neutron flux remains high the capture of neutrons 
may continue until the target nucleus reaches the neutron drip-line. Then the probabil- 
ity of neutron emission increases considerably, and, in the case of U, it dominates over 
the photon emission. The "dynamical" equilibrium between absorption and emission of 
neutrons at the drip-line depends weakly on the neutron density and may last for a long 
time (more than milliseconds) until the /3-decay makes the Z+l nucleus. As the result, 
the drip-line is shifted to the right and absorption of neutrons becomes possible again, 
of course, if the neutron density remains high enough. One can see that the Cameron 
formula leads to similar results for U and Cm nuclei, with saturation at large A. However, 
the MS formula predicts slightly different Cm masses around A 280 as compared to the 
Cameron formula. This leads to a smaller ratio of the neutron and photon evaporation 
widths, r n /T 7 , so that the emission of photons is more probable in the MS case. This 
gives a chance to absorb more neutrons in Cm target, contrary to the U target, and, 
finally, enter the mass region with more bound neutrons. 

Calculation of nuclear abundances: For practical applications, it is very important to 
estimate the amount of new isotopes produced by this mechanism. In Figs. 2 and 3 we 
have presented the evolution of an average mass number of the target nucleus, ignoring the 
distribution around the average. This analysis has provided us with better understanding 
of the main mechanisms responsible for increasing the nuclear mass number under various 
external conditions. As turned out, at temperatures T ~ 10 keV we can simply consider a 
capture of neutrons followed by the emission of photons, without complications induced by 
re-emission of neutrons and fission. In this case we can write the set of coupled equations 
for densities of nuclei with mass numbers A=Aq, Aq+1, Aq+2, ... : 



which can be solved numerically for any given time dependence of the neutron density p n . 
Here v re \ is the neutron- nucleus relative velocity which is close to the neutron velocity v n , 
and a n is the neutron capture cross-section (see section 2), which is assumed here to be 
independent of the nuclide mass numberJl These equations describe the evolution of the 
ensemble of nuclear species due to the neutron capture reactions with gain (A — 1 — > A) 
and loss (A — > A + 1) terms. It is worth noting that final nuclide abundances depend only 
on the time-integrated neutron flux $ = / p n v n dt w v n p n r. According to the analysis of 
refs. [U [10] , typical values of $ reached in low- yield nuclear explosions are in the range 



2 Of course, when the neutron number approaches the drip-line limit, the capture cross section a n will 
strongly decrease, and the solution of the equations will be more complicated. 




(9) 



of 3-10 24 10 25 n/cm 2 . 
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In Fig. 4 we present the distributions of isotope abundances after exposure of a 238 U 
target to a constant neutron flux during fixed time intervals indicated at the lines. The 
figure shows that the distribution of neutron-rich isotopes becomes broader with increasing 
the exposure time, or, equivalently, the integrated neutron flux <3>. According to refs. 
[TO] , in samples of material extracted from underground explosion sites concentrations of 
rare nuclei (A=257) on the level of 10~ n — 10~ 12 were measured by radiochemical methods. 
Even higher sensitivity can be reached by using mass spectroscopic methods (see, e.g., 
[2D]). In typical nuclear explosions mentioned above the observed neutron-rich isotopes 
have not more than 20 absorbed neutrons [7J [10]. According to calculations of Fig. 4, 
where neutron density was fixed at 4-10 22 cm -3 , to produce such isotopes in observable 
amounts we would need r m 0.3/xsec. This corresponds to the integrated neutron flux 
of about 2 • 10 24 n/cm 2 , which within a factor of 2 agrees with the values estimated in 
refs. mHU]. 

From Fig. 4 we conclude that observable concentrations of nuclei with more than 40 
absorbed neutrons can be obtained only if the exposure time is longer than 1/xsec. Ob- 
viously, to increase the number of captured neutrons one should increase the integrated 
neutron flux $, i.e. either by increasing the neutron density p n or by increasing the expo- 
sure time tJH To achieve this goal one can try e.g. to increase the neutron multiplication 
time by introducing neutron absorbers, reflectors, and other construction elements. With 
the exposure time of 3 //sec the fraction of nuclei captured 50 neutrons would be about 
10~ 8 . In this case, even if we take 1 gram of the exotic target material like 248 Cm, the 
concentration of nuclides with A^300 will be 10~ 12 which could be detected by present 
experimental methods. 

These results suggest that the drip-line region can in principle be reached in a prop- 
erly optimized explosive processes. Obviosly, new experiments are required to explore 
this region. They should not only be limited to searching for long-lived superheavy el- 
ements, where the yield may be essentially affected by spontaneous fission. Heavy and 
intermediate-mass nuclei in the vicinity of their neutron drip-lines could be synthesized 
in this way too. For example, we have found that under the same conditions a Pb nu- 
cleus may capture nearly 40 neutrons and come close to the drip-line. This approach to 
the drip-line has an important advantage as compared with heavy-ion reactions, because 
nuclei enter this region with minimal excitation energy. However, due to expected short 
live-times of these nuclei, their properties can only be studied by analizing their decay 
products. 

5. Conclusions 

It is obvious that the nuclear explosions provide the highest fluxes of neutrons under 
terrestrial conditions which will hardly be feasible with any other experimental technique. 
We have demonstrated that irradiation of nuclei with neutrons produced during the first 
microseconds of a nuclear explosion is a very promising way to synthesize neutron-rich 

3 We are taking here about time scales which are still much less than the /3-deacy times, i.e. millisec- 
onds. 
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heavy and super-heavy elements. For comparison, modern spallation sources (e.g. ESS) 
may supply thermal neutrons with densities around 10 12 cm -3 , and nuclear reactors - 
with 10 11 cm -3 [21J. We believe, there exist technical possibilities to increase the time 
and intensity of target exposure in the course of a nuclear explosion. They may include 
preliminary compression of fissile material, introducing neutron reflectors and moderators, 
as well as special construction of targets. The main idea of optimization should be to 
maximize the neutron density and exposure time and to minimize the energy yield of the 
explosion. Another possibility which should be considered in more details is generating 
two or several nuclear explosions within a time delay up to milliseconds in close proximity 
of each other. This delay may be sufficient for very neutron-rich nuclei produced in 
the first explosion to undergo multiple /3-decay. Then the daughter nuclei will be able 
to absorb additional neutrons from the second explosion, and in this way increase the 
resulting mass number, as illustrated in Fig. 1. Even mass production of super- heavy 
elements can be envisaged in the future. 

The methods discussed above can be used not only for synthesizing long-living super- 
heavy elements but also for production of very neutron-rich isotopes in the vicinity of the 
drip-line. These studies are extremely important for understanding the origin of heavy 
elements in the Universe. Hot and dense environments created in nuclear explosions are 
rather similar to conditions in supernova explosions, when heavy nuclei are produced by 
multiple neutron capture in r-process. 

If long-lived superheavy elements are indeed produced in the explosions, an important 
question to be answered is, how to find them among tons of debris. In underground 
explosions the processed material may be extracted and analysed by radio-chemical and 
mass-spectroscopic methods. Generally, the methods used for searching for superheavy 
elements in nature (22] may be applied in this case too. 
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Figure 1: (Colour online) The upper end of the nuclear chart showing the experimentally 
investigated isotopes. The colours stand for the decay modes: yellow - a-decay, red - (3 + or 
electron capture, blue - (3~ decay, green - spontaneous fission. Black boxes are stable nuclei. 
Grey area is the predicted "island of stability". Red arrow represents the heavy-ion fusion 
reactions used so far. Blue arrows show schematically the proposed explosion-based methods to 
reach this island: I indicates fast capture of neutrons during first microseconds of the nuclear 
explosion followed by the /3-decay. II indicates second nuclear explosion within few milliseconds 
after the first one. 
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Figure 2: Evolution of U (top panel) and Cm (bottom panel) mass numbers with time at the 
neutron densities shown at the figure. (The initial target nuclei are 238 U and 248 Cm.) Dotted 
line presents absorption of fast neutrons, other lines - absorption of thermalized neutrons at 
T=10 keV. 
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Figure 3: Effect of different mass formulae on calculations of mass number evolution for 238 U 
and Cm targets. Dashed and solid lines: masses are taken from ref. |16] . dotted and dash- 
dotted lines - ref. |19| . Neutron density, temperature are also shown in the figure. 
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Figure 4: Densities of nuclei produced after capture of neutrons normalized to the initial density 
of 238 U target. Density and temperature of neutron, as well as exposure times (in microseconds) 
are shown in the figure. 
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